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Abstract 

A first-order ordinary differential system with a matrix of periodic coefficients 
Q (y) = Q (y + T) is studied in the context of time-harmonic elastic waves travelling 
with frequency uj in a unidirectionally periodic medium, for which case the monodromy 
matrix M (a;) implies a propagator of the wave field over a period. The main interest 
to the matrix logarithm InM (uj) is owing to the fact that it yields the 'effective' matrix 
Q c ff (uj) of the dynamic-homogenization method. For the typical case of a unimodular 
matrix M (uj) (det M = 1), it is established that the components of InM (uj) diverge as 
(uj — ujq) 1 with uj — > ujq, where ujq is the set of frequencies of the passband/stopband 
crossovers at the edges of the first Brillouin zone. The divergence disappears for a homo- 
geneous medium. Mathematical and physical aspects of this observation are discussed. 
Explicit analytical examples of Q c fj (uj) and of its diverging asymptotics at uj — > ujq 
are provided for a simple model of scalar waves in a two-component periodic structure 
consisting of identical bilayers or layers in spring-mass-spring contact. The case of high 
contrast due to stiff/soft layers or soft springs is elaborated. Special attention in this 
case is given to the asymptotics of Q c fj (uj) near the first stopband that occurs at the 
Brillouin-zone edge at arbitrary low frequency. The link to the quasi-static asymptotics 

of the same Q e g- (uj) near the point uj = is also elucidated. 

Keywords: logarithm of a matrix, ID periodic media, Floquet spectrum, dynamic 
homogenization, high-contrast structure 
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1 Introduction 



The first-order ordinary differential system 

Q(y)ri(y) = ^rj(y) (1) 

with a n x n matrix of continuous or piecewise continuous periodic coefficients Q (y) = 
Q (y + T) is a classical problem arising in miscellaneous models of applied mathematics and 
mathematical physics. Its analysis largely relies on the Floquet theorem asserting that the 
matricant M (y, 0), which is the fundamental solution of ([I]) yielding t] (y) = M (y, 0) t] (0), 
can be factored into the product 

M(y,0) = L(y)exp(iKj/), (2) 

where L (y) = L (y + T) (= M (y n , 0) exp (-iKy n ) with y n = y (mod nT)), L (0) = I (I is 
the identity matrix), and K is a constant matrix [1]. By ([2]), K is defined by the equation 

exp (iKT) = M (T, 0) iKT = In M (T, 0) , (3) 

where M (T, 0) = M is termed the monodromy matrix (its reference to (T, 0) is dropped 
hereafter). It can be calculated by a number of available methods, e.g., using the Peano series 
of multiple integrals of Q(y), or applying polynomial expansion of Q (y) , or discretizing 
Q(y). In many problems the system matrix Q (y) is a continuous function of a certain 
control parameter u, and hence also M = M (uj) . At first glance, the matrix logarithm K (u) 
is well-behaved as long as the logarithm of the eigenvalues q (a;) of M (u) is well-behaved. 
However, it turns out that In M (w) diverges at uj — > ujq, where c^o corresponds to a non- 
semisimple (not diagonalizable) M (c^o) with a degenerate eigenvalue q (uq) whose values, 
being taken on the same Riemann sheet of In q, are situated on the opposite edges of the cut. 
This fairly surprising observation seems to have passed unnoticed in the extensive reference 
literature on the matrix logarithm. The manner in which such a divergence reveals itself 
in the Floquet formalism is discussed in the present paper in the context where Eq. ([1]) is 
associated with time-harmonic elastic waves travelling at frequency to in unidirectionally (ID) 
periodic media. Within this context, the system ([T]) such that consists of n = 2 equations 
and hence is equivalent to Hill's equation describes scalar acoustic (or electromagnetic) waves 
[2113]; the cases where ([1]) consists of n = 4, 6, 8... equations corresponds to coupled waves in 
elastic isotropic or anisotropic media, in piezoelectric or piezomagnetoelectric media, etc. In 
either of these cases, the monodromy matrix M is often called the propagator (of the wave 
field) over the period T. 

The matrix logarithm K (u) (j3J) is a crucial ingredient in the dynamic-homogenization 
approach. Assuming that exp (iK.y) in the ID Floquet theorem ([2]) is a relatively slowly 
varying function, this approach seeks to replace an exact solution M (y, 0) by its 'slow com- 
ponent' exp {iK.y) and hence to replace the actual periodically inhomogeneous material by an 
'homogenized' medium with spatially constant but frequency dispersive properties described 
by the 'effective' matrix 

Q cfr H = tK (a;) , (4) 
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see e.g. [HEJE]. Obviously, the matrix Q eff also provides (regardless of any assumptions) 
an exact solution M (nT, 0) = exp (inKT) at the interfaces between the periods. Another 
aspect of the matrix logarithm K (uj) is related to the Floquet dispersion branches uj (K) or 
K (uj). These are determined by the secular equation for M, 

det [M (uj) - q (uj) I] = 0, (5) 

so that the definition q = e lKT yields %K (uj) T = In q (uj) , or else by the formally equivalent 
secular equation for K, 

det [K (uj) - K (uj) I] = 0. (6) 

The Floquet spectrum is commonly defined over the first Brillouin zone (BZ) Re XT G 
[— 7r, 7r] , which is related to the zeroth Riemann sheet of the single- valued In q = ln\q\+i arg q 
with the cut argg = ±7r corresponding to the BZ edges. The frequency intervals where K is 
real or complex are called passbands and stopbands, respectively. 

The paper is concerned with the typical case where M (uj) is unimodular (det M = 
1) and so the BZ edges contain the passband/stopband crossovers at a set of frequencies 
uj = uo associated with a degenerate pair of eigenvalues q (cuo) of M (cuo) ■ According to the 
background outlined in §2, this is the case for a normal propagation across an arbitrary 
anisotropic periodic structure or for an arbitrary propagation direction in the presence of 
appropriately oriented symmetry plane. The original material of this work consists of two 
parts, § 3 and §4. The first part (§3) deals with the problem in general. It is shown 
that the matrix lnM(o;) = iK (u;) T, and hence Q e g, must have components diverging as 

— 1 /2 

(uj — ujq) 1 when uj — > u , i.e. when the real Floquet branches tend to the BZ edges 
or the complex part of ±K = ir/T + ilmK tends to zero. The eigenspectrum of K (uj) 
certainly remains well-behaved for any uj infinitesimally close to ujq; however, computing the 
Floquet spectrum K (uj) specifically from Eq. ([6]) may become numerically unstable at uj 
close to ujq. A transition is explained from a weakly inhomogeneous to perfectly homogeneous 
elastic medium, for which InM (uj) certainly does not diverge. The second part (§4) presents 
detailed analytical examples of Q c ff (uj) and of its diverging asymptotics for uj — >• uj for the 
shear-horizontal wave in a periodic structure composed of piecewise homogeneous bilayers 
or layers in spring- mass-spring contact. Particular attention is given to the high-contrast 
case with either a soft layer in the bilayer or with a soft spring in the interfacial joint. The 
interest to this case lies in the fact that the first stopband at the BZ edges and hence the 
local divergence of Q c fj (uj) occurs at low frequency that may in principle be made arbitrarily 
small. To this end, a link to the regular asymptotics of the same Q c g (uj) near the point 
uj = is also elucidated. The basic points of the study are summarized in §5. Some technical 
aspects of the derivations of §3 and §4 are detailed in the Appendix. 

2 Background 

Consider elastic waves in a ID-periodic infinite anisotropic non-absorbing medium without 
sources. Choose the periodicity direction as the axis Y and denote the (least) period by T, 
so that the density and the elasticity tensor satisfy p(y) = p(y + T) and c (y) — c (y + T), 
respectively. Take the axis X in the sagittal plane spanned by Y and by the direction to the 
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observation point. Applying Fourier transforms in time and in X brings in the frequency uj 
and wavenumber k x as the (real) parameters of the problem. 

The equation of motion and the linear stress-strain law may be combined into the system 
(CQ) of, generally, six equations. The periodic matrix of coefficients Q (y) defined through 
p (y) , c (y) and uj, k x , is pure imaginary and has the Hamiltonian structure 

Q (y) = TQ T (y) T, (7) 

where the superscript T means transpose and T is the matrix with zero diagonal and identity 
off-diagonal 3x3 blocks (see e.g. [7] for the details). 

In the following we deal with the essentially typical case of a medium with at least a 
single symmetry plane m orthogonal to the axis X or Y. Then the trace of Q (y) is zero 
for any y. Therefore, by the Jacobi identity, M (y, 0) is unimodular and hence so is the 
monodromy matrix M = M (T, 0) , i.e. 

detM = l. (8) 

The identities flTJ) and (JS]) together ensure that for every eigenvalue q a of M, there is a cor- 
responding eigenvalue qp = l/q a where a, (5 — 1, 6. This property has been established in 
[9] for a piecewise constant Q (y) and m _L Y ; its generalization for any piecewise continuous 
Q (y) and for m _L X is obvious. Note that no stipulation of any material symmetry is 
needed if the wave propagates strictly along the periodicity direction Y (i.e. if k x = 0), 
which is when (jSJ) is always true. Also note that the out-of-plane motion with respect to the 
symmetry plane m _L Z of a monoclinic body (which has no other symmetry planes) can be 
cast in the form with property (jSj), see [8]. 

Let uj be a single free dispersion parameter (k x is fixed or expressed through uj). Each 
pair q@ (uj) = l/q a (u) corresponds to a set of dispersion curves K a (uj) = —Kp (uj) in the 
BZ ReK a) pT G [— 7r,7r] , which are symmetric about the line K = 0. In view of ([8]), the 
eigenvalues q = 1 and q = — 1, occurring, respectively, at the centre and edges of the BZ, 
are assuredly degenerate. We are interested in the case q = —1, which is associated with 
a sequence of passband/stopband crossover points at the BZ edges, and specifically in the 
behaviour of the matrix In M = iKT in the vicinity of these points. 



3 Divergence of K (uj) near the BZ edges 
3.1 Derivation 

Denote by uj = ujq the frequency, at which some pair of eigenvalue branches q\ (uj) = l/q 2 (uj) 
of the monodromy matrix M (uj) fall into two-fold degeneracy q\ (ujq) = qi (uJo) = —1 render- 
ing M (ujq) non-semisimple. Consider a function lng = In \q\ + zargg defined on the zeroth 
Riemann sheet with a cut argg = ±7r passing through —1. Let uj lying in the stopband or 
passband tend to ujq from, respectively, above or below. Then qi (uj) and q 2 (uj) tend to e ±l7T , 
thus approaching their degenerate value —1 from the opposite sides of the cut for lng, and, 
correspondingly, lng lj2 (o;) = iKi 2 (uj)T tend to ±Z7T, meaning that two Floquet branches 
tend to the opposite edges of the BZ. 
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This is indeed nothing else than a very standard setup. The state of affairs is, however, not 
so trivial when the same limit u — > ujq is applied to the matrix logarithm In M (uj) = iK. (uj) T. 
It is natural to specify it by asking that both eigenvalues lng 12 (uj) of InM (uj) satisfy the 
above-mentioned definition of In q (the issue of alternative definitions of In M is addressed in 
§3.1 and in §A.2 of Appendix). As we have just observed, these eigenvalues tend to ±in as 
uj — > uj , i.e. they do not approach each other in contrast to the eigenvalues q\ (uj) — > q 2 (uj) 
of M (u;). This signals a singularity of InM (uj) on the path uj — > ujq. 

Let us analyze the local behaviour of InM (uj) for uj = uj + Au; (\Aoj/uq\ < 1). With 
reference to ([8]), denote 

q 12 (uj + Au) « q d ± 5q -> q 12 (uj ) = q d = -1, (9) 

where 6q means the leading-order correction in the small parameter Au/u . For brevity, 
assume the case of 2x2 matrices (the same derivation for the general nxn case is detailed in 
Appendix, § Al). A polynomial formula for a function of a 2x2 matrix M with eigenvalues 
qi 7^ q 2 has a simple form 

/ (M) = g2/ ~ qJ + f fe) ~ f fa) M, (10) 
<?2 - <7i <?2 - Qi 

see e.g. [ID]. Taking §W§ for M (uj + Auj) = M (uj ) + AM with q 1<2 (uj) given by Q yields 

[M(uj ) + AM-q d I} + O(5q), (11) 



/ [M (lu + Aw)] = foi + fo2 I 



foi — f02 ,, , 

V f (qd 



25q 



where /oi,02 = ^ m u;^uj f (Qi,2 (^)) and O is a matrix symbol 'of the order of. For the case 
in hand / = In and /oi,02 = hi (e ±17T ) , whence ( II ip becomes 

(Z7T 1 \ Z7T 

— + - [M (uj ) - q d I] + — AM + O (Sq, AM) . (12) 
Oq qd) Oq 

Since M (ujq) — qdJ- is non-zero for a non-semisimple M (ujq) while 5q tends to zero with 
Auj — > 0, we conclude from Eq. (fl2|) that the matrix logarithm InM (uj) , and thus K (uj) , 
must have components tending to infinity when uj — > co>o- Note in passing that an identically 
zero determinant of the first matrix term on the right-hand side of ( Tl2l) does certainly not 
preclude but, on the contrary, underlies (with due regard for the next term) the necessary 
identity det [InM (uj)] — > it 2 as uj — >• uj . 

Let us now find an asymptotic rate of divergence of InM (uj) in terms of Au; (= uj — uj ) . 

1 /2 

For a non-semisimple M of 2x2 dimension, the leading-order dependence 5q ~ (cj — co>o) 
obviously follows from a quadratic secular equation (jSJ). For the general nxn case, the same 
trend is easy to infer from the leading-order Taylor expansion of D (q, uj) = det [M (uj) — qi] 
about the point of double degeneracy q^ 2 (uj ) = q d , which leads to 

Omitting details (see e.g. [H]), it suffices to note that B is generally non-zero for non- 
semisimple M (ujq) . Thus, by (fl"2l . lnM(w) , and hence K (uj) , diverges as (uj — uj ) ~ 1 ^ 2 
with uj — > ujq. An explicit form of the coefficient B will be exemplified in §4. 
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3.2 Discussion 



A few formal remarks are in order. First it is reiterated that even though the compo- 
nents of the dealt-with matrix lnM(w) = iK. (ou) T diverge as w -> c^o, its eigenvalues 
In gi^ = iKi,2 remain formally well-defined so long as uj ^ uj®. It is also understood 

that the exponential of this In M (a;) at any uj ^ lj certainly reproduces (continuous) M (a;). 
Regarding the infinity of InM (ou) precisely at uj = uj , which is when 5q = on the right- 
hand side of (IT!?!) , it simply tells us that the conventional definition of InM (uj) , which refers 
both eigenvalues In gi )2 (uj) to the zeroth Riemann sheet of In q with the cut arg q = ±7r fixing 
the edges of the BZ lie KT £ [—it, it] , precludes this matrix function of uj from reaching the 
limiting point uj$ of the path uj — > u continuously 

It is clear from the above that shifting the cut in the g-plane away from the point q — — 1 
while keeping In q\^ on the same Riemann sheet leads to a different matrix logarithm In M (uj) 
that has degenerate eigenvalues ln^i (uj ) = lng 2 (^o) an d hence is well-behaved at uj = uj 
and around it. However, this 'gain' for uj near uj is at the expense of one or another 
essential deficiency elsewhere for the redefined In M (uj) . For instance, if the eigenvalues 
lng 12 of InM (uj) are taken on the zeroth Riemann sheet with the cut argg = 0,2tt, then 

— 1/2 

this InM (uj) has the same divergence ~ (uj — uj i) due to the degeneracy gi )2 (ojoi) = 1 
at the set ujqi of passband/stopband crossovers occurring at K = 0, 2ir. An exception is the 
origin point uj = 0, where M = I and so any InM is continuous; however, the low-frequency 
onset of InM defined by taking the cut argg = 0, 2tt has no physical sense (see Appendix, 
§A2). Another possibility is to use a cut argg = <p, <p — 2tt at <p ^ irn, e.g., at ip such that 
< (p < it. Then InM, whose eigenvalues lng 12 = ±iKT lie on the zeroth Riemann sheet, 
is well-behaved with |argg| = \KT\ growing from zero but only until reaching tp, where there 
is a jump to a different matrix InM, for which the eigenvalue lngi has to be shifted from 
arggi = KT to arggi = KT — 2tt with KT > increasing above (p. Note that a similar 
piecewise discontinuity pertains in the BZ Re KT £ [—it, tt] to the logarithm of M that is 
not unimodular (detM ^ 1). Thus, using any 'unconventional' definition of the logarithm 
of M based on shifting the cut from the point q = — 1 is hardly an alternative. 

It remains to settle a natural question concerning the case of a homogeneous elastic 
material, for which the matrix Q is constant, hence M = exp (QT) , and so InM merely 
returns the 'initial' QT, which is certainly continuous in uj. 'Technically', the difference with 
the case of a periodic medium is that a constant Q keeps M (uj ) diagonalizable (semisimple) 
at the degeneracy point qi (uj ) = q 2 (uj ) = —I under discussion^- Assuming M (a> ) = qdl 
in Eq. (|T2|) . its first term turns to zero and thus a continuous InM (ujq + Auj) is defined by 
the second term of (fT2l) . in which AM ~ (uj — uj ) and 5q ~ (uj — uj ) (the latter being due 
to B = in (fT3j) for a semisimple M (ujq) [H]). A transition to (or from) a homogeneous 
material from (or to) a weakly (periodically) inhomogeneous one is also evident: given a 
small parameter e of elastic inhomogeneity, M (u ) — q^I is scaled by e and 5q is scaled by 
(eAa;) 1 ^ 2 , hence, by (fT2|) . the singularity of InM (uj) at uj — > uj$ is proportional to (e/Acu) 1 ^ 2 
and disappears at e = 0. 

In conclusion, let us outline some exceptional cases that are theoretically possible due to 

1 For a constant Q, this degeneracy of gi,2 = e ikyT implies nothing more than an odd number of half- 
wavelengths within the interval Ay = T - note no relevance to degenerate eigenvalues k y of Q that do render 
Q and hence M = exp (QT) non-semisimple. 
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'incidental' occurrence of M (u; ) in a peculiar form. First, a non-semisimple M (u ) does not 
preclude vanishing of the leading-order coefficient B ffT3"l) 9 [TTj ; if it happens to be zero then 
{Sq) 2 is given by the higher-order terms of the Taylor series of D (q,ui) about c^ch i n which 
case Eq. (fT2"|) (where M (uj ) ^ q d I) leads to In M (lu) ~ (lu — lu ) ~ m l 2 with an integer m > 2. 
Secondly, a (periodically) inhomogeneous medium M (lu ) does not rule out a possibility for 
M (lu ) at a degeneracy point to remain semisimple (such an option is usually associated 
with a stopband of zero width). Finally, a semisimple M (luq) may, in principle, also cause 
diverging InM (lu + Alu) - it is the case when 5q ~ (lu — u; ) m with m > due to 
incidentally vanishing higher-order derivatives d 2 D/du> 2 , d 2 D / dqdu etc. in the Taylor series 
of D (q, lo) about luq, whence In M (ou) for lu — > luq diverges owing to the term (Sq)~ x AM ~ 
(^-^ )" m/2 in Eq. (TTJ]). 



4 Examples of Q e g = zK 
4.1 Bilayered unit cell 

This section is intended to illuminate the preceding general development by way of its ap- 
plication to simple examples of a scalar acoustic wave in a periodically repeated sequence 
of pairs of homogeneous layers. With this purpose, we first remind the 2x2 setup for an 
arbitrary ID-periodic medium [2], [3] and detail the formulas describing the 'effective' matrix 
Q c fj = iK for this framework. Then we further elaborate Q e ff for the case of a bilayered unit 
cell. 

4.1.1 2x2 setup 

Consider a 2x2 unimodular monodromy matrix M (u) . Its eigenvalues 

q l 2 = ItraceM ± R, where R = - J (trace M) 2 - 4 ( = qi ~ q A ; (14) 
define the Floquet wavenumbers 

zA'i^T = ±iKT = In = izarccos ^-traceM^j = ±2zarccos ^- Vtrace M + ; (15) 
and the equation 

traceM = -2 (16) 

defines the set of frequencies uo = Uo of passband/stopband crossovers at the BZ edges 
KT = ±7T where q x (u ) = q 2 (uj ) = q d = —1, see |2] [3]. 

Introduce the 2x2 'effective' matrix Q c g = zK, which is related to M by the equality 
M = exp (zKT) and which has eigenvalues (fl~5l) understood under the standard definition of 
the functions In and arccos, so that Re KT G [— 7r, 7r] . Then Eq. fTTOl) specified for / (M) 
1 n M gives 

(17) 



^ = ~R 



M-- (trace M) I 
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The same result may certainly be obtained by equating M to exp (z'KT) , which follows from 
the same (ITUi) (re- adjusted to / (K)) in the form 



exp(iKT) = (cos KT) I 



sin KT \ \, i? 
i — - — ) K =- (trace M) H K 



K 



K 



due to using the condition K\^ = ±K equivalent to fixing the appropriate definition of 
matrix logarithm InM. 

Consider now a vicinity of the BZ edge. Eqs. ( fT4"l) . (1T51) expand in small Au = uj — ujq as 



9l,2 \U) 



where it is denoted 



■1 ± VBAuj + O (Au) , K(u)T 



ix + iy/BAu + O (Au) 



B 



BAu + O [(Aw) 2 ]. 



d 

- — trace JVl 
do; 



(19) 



(20) 



which is non-zero for a non-semisimple M (o>o) (barring the theoretical exceptions mentioned 
in the end of §3.2). Inserting ([19]) and in (fTTjl yields 



/7T 



v^BAU + O (Au) 



A 



dM 

da; 



Ira 

2 



Au; + O (Aa;) 



(21) 



'fiAa; + O [(Aa;) 2 ] 

where A denotes a non-zero nilpotent matrix 

A = M(u; )-g d I = M(u;o)+I (A 2 = 0). (22) 

Eq. d2IJ elaborates (El (with due regard for AM /Sq ~ O (5g)). Note also that Eq. (jT5]| q for 
i?, defined in ( |T4l) 9. to leading order reads 5q = y/BAu which is recognized as the equation 
(|T3|) i . Correspondingly, the definition (|20|) of the coefficient B is equivalent to Eq. (fT3l) 9. 
which specializes for the given case (of 2x2 M with qa = — 1 at o; ) as 



5 = - 



do; 



det (M - ql) 



trace 



A 



dM 

do; 



(23) 



Expansion (121!) shows that the 'effective' matrix Q c g (a;) has well-behaved eigenvalues 
±iK (u) — > ±iir/T at uj — > Uq, while its components diverge due to non-zero A with a 



common factor 



{u - u; J 



-1/2 



. It is also seen from Eqs. (pT]) - (T23|) that A and B for a 
weakly inhomogeneous unit cell can in general be scaled by the same small parameter e (= 

1 /2 

for a homogeneous limit), and so the singularity of Q e g (uj) at cu — > u is scaled by (e/ Au) 1 
as argued in §3. 

4.1.2 Q cff for a bilayered unit cell 

Let us narrow our analysis to the case of a two-component piecewise constant unit cell. 
Specifically, we consider the shear horizontal (SH) wave in a periodic structure of perfectly 
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bonded pairs of isotropic homogeneous infinite layers j = 1,2, each with constant density pj, 
shear modulus /x - and thickness dj. For the sake of the brevity of explicit formulas, assume 
the wave u (y) propagating along the axis Y normal to the interfaces (k x = 0). Hooke's law 
a (y) = fjbju' (y) and the equation of motion a' (y) = —pjU 2 u (y) combine into the system 
(CQ) with the state vector r\ (y) = (iouu, cr) T and the piecewise-constant periodic 2x2 system 
matrix 

Q^iusJ ^ Z 'i ) y j = 1,2, (24) 

which leads to the propagator M (T, 0) = e Ci2d2 e (:ildl = M (u) through the period T = d\ + d 2 
(the monodromy matrix) in the form 

, , . . / cos ip? cos ib* — f 1 sin ib~ sin xb* cos-0 9 sin 7/^ + sin-0 9 cosipi \ . , 

M(u)=[ .„ , . , z * , Zl . 1 , z f 2 . , . . 1 , (25) 

\_ %Z\ cos ip 2 sin ip 1 + iZ 2 sin ip 2 cos y i cos^ 2 cosyj 1 — -rr sin^ 2 sm-^ 1 y 



where Sj = yPj/Pj is the slowness, Zj = *Jp~Jp~j the impedance and ^ = wSjdj the phase 
shift over a layer. Passing in (125]) to an oblique propagation amounts to merely premultiplying 



ibj and Zj by yl — with a fixed s x = k x /u>. Inserting M into the basic relations (1141 - 

(j!6p provides the textbook equations implicitly defining the Floquet spectrum u (K) and its 
stopband bounds to = u at the BG edge for a bilayered unit cell, e.g. [2]. 

The 2x2 'effective' matrix Q c fr = ?K for a bilayered unit cell follows from (Ti~7j) and ( |25l) 
in the form 



Qeflf (w) 



z'if / - \ ~ J|) s in ^2 sin ^1 cos 2 sin ^ + ^ sin ^ 2 cos ib x 



R yiZ\ cos "02 sin ibi ~\~ iZ 2 sin ip 2 cos "0^ | — J 2 -^ sin-0 2 sin-0 1 

(26) 

It is easy to check that the eigenvalues of this matrix are ±iK, and that it reduces to 
( H/^i , when Zi = Z 2 , Si = s 2 . As another consistency test, we note that (I2"B"1) provides the 
well-known low-frequency asymptotics of Q e ff, whose diagonal and off-diagonal components 
expand in, respectively, even and odd powers of iu as follows: 

Qeff {U) U/U0<X = (Q) + ^ (Q2Q1 " Q1Q2) + 

where 

(Q (c)) = Q,f + Q 2 f , (^-i) = + JL* 

(28) 



\P) - Pit + ft T ' K - ^2- ^— Mi 



Additional explicit insight is gained by noticing that trace M + 2 with M given by 
can be factored as 

trace M + 2 = /+/_, 
U = vfe + Z *) cos ^ ± - 3,) cos ^ 



(29) 
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whence Eqs. ( ED , (jl5]) provide 



itf (w) T = In ( - 1 + r) = 2% arccos ^ , # = *//+/_ ( J ^ L - - 1 ) , 




and Eq. f|T6|) takes the form 

/+/- = 0, (31) 

showing that the set uj = ujq consists of two families given by zeros of f±. Evidently, this split 
reveals the symmetric/antisymmetric decoupling of the problem. As a result, the expansion 
fl2T]) of Q cff (u) about the points uj = lo , when applied to the matrix Q cff fl26|) in hand, 
admits compact formulas for its leading-order parameters B (120]) and A (122]) as follows: 

* = - {f^) U0 = (f - f ) sin ^ + ^ 2 Bin^) , 

COS V'l + COS ^ 2 ^ S ^ 1 ' 1 i J l ~ ^ S ^ rl i J 2 \ 

iZ 2 sin ?/> 2 — iZ\ sin , 1 — cos — cos ?/> 2 



A = ± 



where ^ - = uj Sjdj are referred to cl> , and the upper or lower sign corresponds to / + = or 
/_ = in (l3~il . respectively (see §A3 of Appendix for derivation of (l32]) o) . The derivative 
(dM/da>) , which also appears in ( |2T]) . can be obtained due to M = e^^e^ 1 ^ in the form 
expressed through the matrices Qj ( l24j) and A as 

^) = — (d 2 Q 2 M + d 1 MQ 1 ) ulQ = — [^Q 2 (w ) A + ^AQ X ( Wfl ) - (Q (w )> 
duj J £Jq wo L T T 

(33) 

Its plugging in f|23]) o and taking note of A 2 = yields another definition of the coefficient 
B, 

T 

B = trace [A (Q (co ))) , (34) 

UJ 

which for the given case of a bilayered unit cell is equivalent to (120]) and (123]) . It is easy to 
verify that fljgj) with (PJHD and (J32D 2 leads to (j52jli. 

The following analysis for highly contrasting layers and for layers in spring-mass-spring 
contact makes an extensive use of the factorization (129]) and the consequent formulas. 



4.1.3 High-contrast case 

It is instructive to specialize the above considerations to the case of high contrast between 
the material properties of two layers composing the unit cell. Suppose that, e.g., the second 
layer is much softer than the first one: 

yU 2 /yU 1 = e 2 (=> s 2 ~ Z 2 ~ e) , where 0<e<l. (35) 

The main interest of the high-contrast case is that the first stopband at the BZ edge occurs 
in the low-frequency range, which is scaled by e and implies ip ± = O (e) , ip 2 = O (1) • I n this 
range, the propagator (125]) is approximated to leading order in e as 

■\A( A -( cos^ 2 -/5sin^ 2 ^sin^ 2 \ ^ 

{ U) ^o(e)- \ iz 2 (simj; 2 + f3 cos *j; 2 ) cos^ J 1 1 ' 
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where 



Zo 



UJ- 



and Eq. (fT6|) with M (]36j) defines the stopband bounds to = uoq by 



cos^ 2 - -sm^ 2 



^2 £ • "0S 



1 -v^ cos — cos sin — 



0. 



(37) 



(3* 



The latter, factorized, form is Eq. (I3T1) with approximate f± ( 129]) -?. So the first stopband 
is bounded by the least roots of /+ = and /_ = which, to leading order in e, are the 
first zeros of the cofactors of (158]) ■>. The upper bound corresponding to /+ = is close to 
the first thickness resonance ip 2 ( = ^2^2) = ^ °f the s °ft layer. Denote the lower bound 
corresponding to /_ = by Q (= minc^o) . It is approximated by the least root of equation 



tan(^ 2 /2) = 2//3, 



(39) 



which involves coupling of the layers. Note in passing resemblance and dissimilarity between 
this simple model (see also §4.2) and the textbook case of a high-contrast diatomic lattice 

With a view to highlight the low-frequency behaviour of Q c g (u;), let us focus our attention 
on u) ranging from uo rs Q and going down the first Floquet branch to uo = 0. Substituting 
(USD in d2SD yields 



iK 



■f sin ^ 2 



iZ 2 (sin ip 2 + P cos '02; 



Qcff M 

where by §U§ and ([38]) 

iK (ou) T = In (cos ^ 2 — f s i n ^2 + ^) = 2z arccos 



cos 



^sin^ 2 
f sin ip 2 



cos ^ — I sin ^ 



(40) 



R(u) 



2 sin ( f cos 9 + s i n 9 



f sinf 



COS -y- 



1/2 



1/2 



(41) 



The singular term for Q e g (o>) (]4"0]) as u tends to the first stopband bound Q is Q fr (to) oc 
^7== A (see (J2U) with yf^j = i^VL-u and 



B 



1 



[ / g(V' 2 + sin^ 2 )] 6 



A 



1 2j 

1 Z 2 (3 



(2/(3f \ *Z 2 p/2 



(42) 



Eq. (142D follows from ( 132]) which is taken with the lower sign (since f2 is defined by /_ =0) 
and confined to leading order in e (in accordance with the accuracy of (136]) and hence of 
(l4"0"]) ). The asymptotics of the same Q c ff (u) (I4"U1) near the origin point u; = is given by Eq. 
(127]) with 

/ -1\ 1 ^2 Pi / j q\ 

(j. 2 T T 2 fi 2 

which also implies taking leading order in the high-contrast parameter e. Note that B pro- 
vided in (|32Di satisfies Eq. (PJ with (Q) given by ([43]). 
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4.2 Layers in spring- mass-spring contact 

As another example, consider propagation of the SH wave through a structure of identical 
layers of thickness T in spring- mass-spring contact. Denote the rigidity of each of two 
springs by 7 and the mass by m. Note that the physical dimension of m is voluminal density 
times length. The monodromy matrix M = M (T, 0) for the state vector r] = (iurn, a) T is 



M = M int M;, where M; = exp (QT) is the propagator across the layer with Q given by 
(no subscripts j = 1,2), and M int is the propagator across the spring-mass- spring interface: 



y uj 2 rn 2iu> f -y ui 2 m \ \ / -y 2uj 2 2iu) i -y 



..2 



*j \ v -^ 7 n 4 '.-r h (44) 

where Q r = ^/27/m is the resonant frequency of this joint. Thus 



MM 



1 - %- cosV - ^ 1 - & sin^ i 1 - %- sin^ + ^ 1 - % ) cos^ 



(l - f£) sin^ + ^ cosV ; (1 _ *fj cos ^ _ |g sin ^ 

(45) 

A factorized form fl3TT) of the equation ffT6]) defining the stopbands at the edge of the BZ 
holds with 

/+ = 4 [( 1 -^ 2 ^) cosf-^sim 



/- = cos I - ^sinf, 



(46) 



where u 2 /Vl 2 = i/j 2 mp/2p r yT 2 is used to write f± as functions of the phase shift tjj = ujT^J ~p[p. 
It is seen that / + (ip) depends on both spring and mass parameters 7T//X and m/pT, while 
/_ (ip) depends on p/'jT only. 

Let us again specialize our consideration to the high-contrast case of a similar 'stiff/soft' 
nature, now by assuming a relatively small rigidity 

jT/p < 1 (47) 

of the springs supporting the mass. Like before, we are interested in the first stopband at the 
BZ edge. Given ( )47|) . the least roots tp ± of f± — and the corresponding stopband bounds 
Q± to leading order are 



^+ = min (v / ¥f > n ) =► Q + = min ^ = min f 



(4* 



1" 



where f2; is the frequency of the thickness resonance of the layer. The question is which of Q + 
and Q- is the lower frequency bound. Since Vl 2 _/VL 2 = m/pT, it is evident that a heavy mass 
m ^> pT ensures Q + = Q r < fi_; a 'medium heavy' mass m ~ pT implies commensurate 
Q + = Q r ~ fi_, and a light mass m <C pT ensures < For the two former cases, 
the whole first stopband is confined to the low-frequency range in the sense that both its 
bounds provide a small phase ij) <C 1. In the latter case of a light mass, decreasing the small 
parameter m/pT keeps the lower bound at ip_ <C 1 and lifts the upper bound up until the 
phase ip + reaches n, i.e. Q r reaches 
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Consider the range 1/) < 1 containing one or both bounds or fi_ and Q + = Q r , 
respectively) of the first stopband at the BZ edge. Expanding (1451) to leading order in small 
if), bearing in mind (14 7p . and using the notations (148 p of Q_ and Q r yields 



A ^Uf i 

which observes det M = 1. Inserting (149]) in (TTTj) gives 
in which 




i-fT (a;) T = In (2a — 1 + i?) = 2i arccos \fot, R (uj) = 2i\J a (1 — a) x 
with a = (1 - cu 2 /ft 2 ) (1 - u 2 /nl ) . [b } 

The latter follows from a similar expansion of f± (if)) (146D at ^ <C 1 and 7T/ yU <C 1, approx- 
imating the l.h.s. of Eq. (|3T|) as /+/_ = 4a and plugging it into (130]) . 

According to (150]) and fl5~T]) . the matrix Q cff (uj) , as expected, experiences the square-root 
singularity at the BZ edge; however, it does so in a different way when uj approaches either 
f2_ (light mass) or Q r (if Q r fulfils if) -C 1; heavy mass). By (j4"9"|) and (150]) . all components 
of the matrix A = M — I are non-zero and hence all components of Q e ff (uj) diverge 
when uj — > Q_ . This is a typical option for a singularity of Q e ff (a>) . On the other hand, 
A = M (Q r ) — I has only left off-diagonal component being non-zero, and hence only this 
component of Q c g (uj) diverges when uj — » Q r while the others tend to zero. This is rather 
an unusual option, which is due to the approximations underlying a simple form ([49]) and 
( I50p of M and Q c fr- The transition between the two above options occurs at Q_ = Q r (i.e. 
m = pT), in which case uj$ = fi_ = Q r implies the stopband of zero width that yields 
a semisimple M (uj ) = —I so that Q c fj (uj) is well-behaved at uj — > uj (it is one of the 
extraordinary possibilities mentioned in the end of §3.2). For either of these cases, the low- 
frequency asymptotics of Q c fj (uj) (150]) is given by (T27]) with the effective properties taken to 
leading order in the soft-spring parameter 047p . i.e., with 

(A* l ) = ^ (p)=P+ ¥ , «=^, (52) 



where 7 S = 7/2 is the rigidity of two identical springs in series (cf. (T43]) ). 

The two types of singular behaviour of the 'effective' matrix Q c fj (uj) defined by (130]) . 
(I5ip are illustrated in Fig. 1. It displays the off-diagonal components normalized by their 
statically-averaged values (Q)^ (~ uj, see (12"7]) with (|52]) i 9) and compares the diagonal 
components to their leading low-frequency term (~ uj 2 , see ( ]27]) with (]52]U ). Specifically, 
the plotted curves are defined as (x) = (Qcfr)^- («j = 12, 21) and y i4 (x) = T (Qcfr)^ 

(2/22 = —yn) with x = uj/VL- when fi 2 , = (1/3) f2 2 (Fig. la) and with x = uj/Q r when 
tt 2 _ = 3fi 2 (Fig. lb), where fi 2 /fi 2 = m/pT (> 7T/2/1) . 
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0.5 1 0.5 1 



(a) (b) 

Figure 1: Frequency dependence of components of the 'effective' matrix Q c fj (1501) in the first 
passband, which is (a) uj G [0,fi_] with Q'i/Q 2 . = 1/3 and (b) uj G [0, Q, r ] with fi^/f^ = 3. 
Black curves are the off-diagonal components ij = 12, 21 normalized by their statically av- 
eraged values; grey curves are diagonal components, whose leading low-frequency evaluation 
(~ uj 2 ) is shown by dashed line. The curves definition is specified in the text. 

Note in conclusion that passing to the case of an oblique propagation (k x = ujs x ^ 0, see 
note to (T23|) ) implies replacing the entries of layer density p by p — s 2 p. Moreover, this case 
enables further 'ramification' of the spring-mass-spring model by means of recasting the point 
mass m as an 'elastic' mass m (1 — c^s 2 ) with its own shear velocity ct (then VL 2 becomes 
Q 2 = 27/m (1 — c^sl)). It is also noted that the case of layers in 'pure spring' contact (i.e. 
without a mass) is described by the above formulas taken with m = (Q 2 — > 00) and with 
7 S = 7/2 as the rigidity of the spring joint, or else by the formulas of §4.1 taken in the limit 
d 2 — > 0, /i 2 — > while keeping 7 S = p 2 /^2 finite. 

5 Summary 

Components of the matrix logarithm In M, where M = M (T, 0) is a unimodular propa- 
gator matrix relating the acoustic wave field with a frequency uj at one and the other 
ends of a period T of ID-periodic anisotropic medium, have been shown to diverge when 
the frequency u tends to the values c^o of passband/stopband crossovers occurring at the 
edge of the first Brillouin zone (BZ). Explicit analytical examples of the 'effective' matrix 
Q cff (uj) (= iK(u)) = |;lnM(a;) and of its diverging asymptotics near the BZ edges were 
provided for the simple case of a scalar waves in a two-component periodic structure of sev- 
eral types, including its high-contrast model when the least of uiq may be made arbitrarily 
small. 

Whereas the components of matrix Q e fj diverge at u) — > Uq, it is understood that Q c g 
for any uj 7^ Uo yields a continuous M = exp (Q e gT) and has a continuous eigenspectrum 
which is in one-to-one correspondence with that of M. Thus, invoking a diverging Q c g 
for formulating a time-harmonic wave propagation through a finite or infinite number of 
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periods cannot create any difficulty, because this phenomenon can be fully described via M 
and its eigenspectrum. At the same time, divergence of components of Q e fj calls for careful 
interpretation if the governing system (pQ) is taken with Q c g in place of the actual matrix 
of coefficients Q (y) and is then viewed in the same sense as the 'true' system ([T]), i.e., as 
incorporating the equation of motion and the constitutive law, but now with the constant 
coefficients Q c g (u) of the fictitious homogenized medium. 

Explicit results of this paper can readily be adjusted to other physical problems whose 
mathematical formulation admits reduction to Eq. ([T]), see e.g. [12] • Further development is 
underway to analyze the Floquet dispersion near an arbitrary passband/stopband crossover 
occurring anywhere in the BZ of ID-periodic structure. Another interest lies in the potential 
extension of the analytical means of the paper to more complicated cases, like in [13] , whose 
exact mathematical statement does not reduce to jl]). 
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APPENDIX 

Al. On the divergence of the logarithm of a n x n matrix M(cu) 

Let M(cj) be a n x n non-singular (det M^0) matrix, continuous in uj, with eigenvalues 
qj{uj). Denote M(cj ) = M , qj(uj Q ) = q® and suppose that q® = q 2 = qd while all other g° 
(j ^ 1, 2) are distinct. Consider a small neighbourhood of uj , where 

M(uj) = M + o(l), Qj (u) = q" + o(l), (53) 

and all qj(ui) (j = 1,2, ...,n) are distinct. Assume that the matrix M with a degenerate 
eigenvalue qd is non-semisimple, i.e., that the Jordan form J of M is 

J = P©S with P = g d I 2 + R, R= f ° J V S = diag(g 3 °,...,g°), (54) 

where I m denotes the m x m identity matrix. Thus the spectral decomposition of M(u;) is 

M(w) = C (w) diag (q 1 (uj) , q n (u)) C" 1 (w) for u ^ co , . , 

M = C J Co 1 for = u , { 1 

where C (u) and Co are matrices whose columns are linear independent eigenvectors of, 
respectively, M(c<j) and Mo (note that Co, which includes a generalized eigenvector of Mo, 
is certainly not C (uo) , which is singular). 

Introduce a logarithm of M (uj) with uj ^ ujq, 

lnM(cj) = C (uj) diag (In q 1 (u) , ...,\nq n (uj)) C' 1 (uj) , , . 

where lng^ = In \ qj \ + i(argqj + 2irkj), kj £ Z. 

This is a general definition in the sense that, while observing indeed the equality exp [In M (uj)] 
M (uj) , it permits taking each Inqj in ( l56l> ) on any fcj-th Riemann sheet. Let us further sup- 
pose that 

IngO-lngO = 2tu, (57) 

which implies either that gi (a;) and q2 (uj) tending to qd as uj — > uj are defined on adjacent 
Riemann sheets (k\ — k 2 = 1) and is away from the cut, or, alternatively, that qi (uj) and 
q 2 (uj) are taken on the same Riemann sheet (ki = k 2 in ( I5B1> )) with the cut such that q® and 
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q® are located on its opposite edges. The latter option with /c 1)2 = is directly related to 
the physical context discussed in this paper. 

Our purpose is to show that, under the aforementioned assumptions, the asymptotics of 
In M (uj) at uj — > ujq is 



In M (uj) 



2ixi 



rA + O 



q 1 (Lu)-q 2 (uj) Wl(^)-92(w) 

A = C [J - diag (q d , q d , g 3 °, q° n )} = C (R 



with 

0n-2)C Q 1 7^ n , 



(5* 



where m ismxm zero matrix and the other entries have been defined above. 

The derivation of (1581) is based on the Lagrange-Sylvester formula [10] with due regard 
for (I53p . (1551 ) and (|57|) . Along these lines, we manipulate InM (ou) as follows (omitting for 
brevity the argument uj of M (uj) and qj (uj)): 



ELi n 



M -9?I 



o 



n 



i>3 

2ni 
91-92 



M0-9°In 

9d-9? 



In g£ + o 



9fe~9j 

^"^" (lngg -lng° 



91-92 
1 

91-92 

+ o 



M - 9j I r 



2?n 

91-92 



Cr 



ix> 3 
rij>3 



9d~9 



(M - q d I 

Jo _ qdln) 



o 



91-92 
1 



91-92 



Co" 1 



o 



1 



91-92 



Next we invoke ( 1541) and observe that 



(59) 



n 

\j>3 



q^n 



9? 



R©0 



n-2j 



(60) 



which is due to 



-9^2 



>3 



9d-9 



I (P - q d h) 



n 



(i 2 + e;= 

S-gPl»-2 
i>3 9d-9? 



n 



9d-9, 



i>3 

r) R 



9d-9, 



rR 



(S - gj n _ 



R; 

= 0„_ 2 . 



R 



(61) 



Note that an essential simplification of (1601) is a consequence of R 2 = O2, yielding (IBTli ). 
Finally, inserting (16"0~|) into (I59p delivers the sought result (1581) . Admitting A = n in 



would lead to a contradiction 0„ = Cq 1 ACo = Ro © n _ 2 7^ n , hence A 7^ n . ■ 

Equation (|58|) shows that the condition (1571) leads to divergence of In M (uj) with qi (uj) — ► 
g 2 (w) at a; — > c^o- For a unimodular M, taking (1581) with gi (u;)— g 2 (w) ~ 25g gives InM (uj) = 

|^A + o • in the case of 2x2 matrices, A = CoRCq 1 = Mo — q d I and hence (158|) provides 

the leading-order term on the right-hand side of (112)) . 

A2. Low-frequency asymptotics of InM defined over the Brillouin zone [0, 2ir] 

Interest in the 'effective' matrix Q e fj = zK = ^ In M is often confined to the frequency 
range uj G [0, Q] occupied by the first passband, i.e. by the first Floquet branch. The 
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logarithm of a unimodular M does not diverge at uj — > f2 if, contrary to the conventional 
definition, its eigenvalues lng are defined on the zeroth Riemann sheet with a cut argg = 
0, 27T. Like any other InM, it is also continuous for uj — > 0. We will, however, demonstrate 
that its low-frequency asymptotics has no physical sense and thus the so defined In M (uj) is 
of little if any practical value. 

For brevity, consider the case of a 2x2 matrix Q (y) given by (T24T) . in which, however, we 
keep arbitrary periodic p(y) , p (y) instead of p •, Py The matrix M (T, 0) = M expands as 
the power series 



I + iujT 



M = I + f T Q (y) dy + J™ Q (y) Q ( yi ) dyd Vl + ... 

(p) o i + 2 ^ Mo ( P ) (p- 1 ) 



+ ..., 



(62) 



where (■) = (•) d? and k = f £ [p (?) p" 1 (?i) - p' 1 (?) p (?i)] d?d?i. For an oblique prop- 
agation {k x = ujs x 7^ 0), p should be pre-multiplied by 1 — s x p/ p. If the period T consists of 
two homogeneous layers, then (p) , (p" 1 ) and k reduce to (|28|) . 

Reserving the notation InM for the conventionally defined logarithm of M, introduce 
another logarithm InM with the aforementioned 'modified' definition, so that 



InM 
InM 



Cdiag (lngi, lng 2 ) C 1 with lng 
Cdiag ( ln<7i, lng 2 ) C _1 with lng 



In \q\ + i argg, — ir < argg < 7r; 
In \q\ + i argg, < argg < 2tt, 



(63) 



where gi.2 are eigenvalues of M, and C is a matrix of eigenvectors of M. Obviously, taking 
exp of both In M and InM returns M. However, these two matrix logarithms are essentially 
different. Note that the standard definition used in ( l63l) i allows the Taylor series ln(l + z) = 



|z 2 + ... for z < 1, whereas In (1 + z) used in ( 1631) 9 is not analytical near z = and 
hence does not admit the Taylor expansion. This underlies a drastic disparity between the 
low-frequency asymptotics of InM and InM. 

For small uj, when g 12 {uj) are close to 1, InM and InM are related as follows 



InM = Cdiag In q u \nq 2 C 



In M + Cdiag (0, 2ni) C" 1 = In M + -g^ (M 



Cdiag (In qx, In q 2 + 2ni) C 1 



(64) 



where, with reference to (TT4T) and (162 



g li2 (uj) = 1 ± iuT^/{p) (p-i) + l - (iujT) 2 (p) (p- 1 ) + O {uj 3 ) . 



(65) 



Hence an explicit difference between InM and In M at uj — > is 



2m 



2-iri 



iujT 




92-91 ~~ 2iLuT v /(p)(n- 1 )+0(w 2 ) 



0(^ 2 



Ofcu 2 



(66) 
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A 



The low-frequency asymptotics of In M = Q e fjT readily follows from (162]) on the basis of 
the Taylor series of ln(l + z), see its example (12 7|) . It has a perfectly clear physical meaning, 
for Q e fj tends to zero when uo — > 0, and to an appropriate matrix Q of a homogeneous 
medium when the inhomogeneity tends to zero (cf. (1271) and fpMl) ). As regards InM, Eqs. 
flM|) and (1661) show that its discrepancy with InM is non-zero even at uj = 0. Thus, contrary 
to InM, the asymptotics of InM near u = has no physical sense. 
A3. Explicit form (132 j) 9 of the matrix A and its properties 

Consider the matrix A = M (lj ) — q d I = M (u ) + I defined at the BZ edge, see (T2"2~]) . 
Substituting the propagator M through a bilayered unit cell given by (125]) leads to 

1 / — | [Z\ — Z 2 ) smtp 1 sin ^2 i (Z 2 smip 1 cos^ 2 + Z\ smip 2 cos^ 1 ) 

Z V Z 2 V iZ]_Z 2 (Z\ smijj 1 cosi/j 2 + Z 2 sin?/; 2 cos^) | (Zf — Zf) sin-?/^ sin^ 2 

(67) 

where ip^ = ojQSjdj (j = 1,2) and u is implicitly determined by Eq. (|T6|) or its equivalent 
( ]3T]) . In the following, the reference to uo = uj will be understood. The objective is to 
manipulate (16 7p into a form that is transparent. 
Introduce the auxiliary notations 

Z± = Z 1 ± Z 2 , ^± = ^(^i±^ 2 ), a ± = Z ± cos^ ± , b± = ^(sin^ + ± sin^_). (68) 

Note the trigonometric identities 

sin ip l sin tp 2 = 4&+&_ = cos 2 — cos 2 
Z 2 sin tp l cos ip 2 + Z\ sin ^ 2 cos tp 1 = a + sin ip + — a_ sin ^_ ; (69) 
Z\ sin cos ip 2 + Z 2 sin ip 2 cos ^ = a + sin ip + + a_ sin ip_ . 

Next we use Eq. (13 lj) . which defines two families of the stopband bounds c^o given by either 
/ + = or /_ = 0, i.e. by either a + = —a_ or a + = a_ (see ( 1291) 9 and (165]) ). Combining 
these equations with (16111) leads to the following alternative expressions for the diagonal and 
off-diagonal elements of A: 

Z Z 2 

A n = -A 22 = -2 + — b+b^ = ±——a + a^ = ±(cos^i + cos^ 2 ), (70) 
%2 il 

Au = -7T-^a±b± = T-^-a T b±, A 21 = i2a±b T = =p2a T 6 T , (71) 



and 



where 



2a T b T = Z\ sin ip-^ — Z 2 sin ip 2 , 2a±b± = ± (Z 2 sin ip 1 — Z\ sin ip 2 ) . (72) 

Except for the first expression in ( 1701) . all others may be called conditional as they depend on 
which of the two families of c^o they are referred to. The compact form of these expressions in 
(17D]) -(172 ]) implies that the upper/lower signs and, simultaneously, the upper/lower subscripts 
are related to / + = and to /_ = 0, respectively. By using these expressions, Eq. (16Tj) can 
be recast in the form 



A = ±( cos^i + cos^a ^sin^i - ^sin^ 2 



1Z2 sin ip 2 — %Z\ sin if) 1 — cos tp l — cosip 2 
± (e Qad2 G + Ge Qia!l ) , where G = ( J _° x j 



(73) 
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with ± corresponding to f± = as above. This is Eq. (1521) 9 presented in §4.1.2. 

By the definition, A = M (u> ) — qdX = for a homogeneous medium, in which case (1751) 
holds with Zi = Z 2 , Si = s 2 and with ip 1 + ip 2 = tt (2n + 1) due to w = u . The matrix 
A for a periodically bilayered medium may incidentally vanish if both cos and cos ip_ at 
uj = ujq happen to turn to zero at once, i.e., if ip 1 and ip 2 i R dZ3j differ by ±7r and in addition 
one of 2 is equal to 2im. In general, A is non-semisimple with a zero eigenvalue and hence 
it must also admit a dyadic representation via its null vector u. This representation further 
specifies due to the identity M _1 = TM + T following from ([7j), which may also be combined 
with the material-symmetry relation M = GM*G to give M- 1 = JM T J- 1 , where J = TG; 
T is a matrix with zero diagonal and unit off-diagonal elements; * means complex conjugate 
and + Hermitian adjoint. Hence 



A = u ® v (= UiVj) , where Au = 0, v = Ju = Tu*, nfli = 0. 



(74) 



Indeed, Eq. (1731) may be re-arranged in the form 



A 



iZi sin ip 1 - iZ 2 sin ip 2 \iZ 2 sin tj; 2 - %Z X sin 



±1 / cos^i + cosip 2 



) 



( 



iZ\ sin — iZ 2 sin ip 2 

cos^/'i + cost(j 2 



) 



(75) 



which satisfies ( 1741) . 
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